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ABSTRACT 

We have proposed a method to calculate the relative magnetic helicity in a finite 
volume as given the magnetic field in the former paper (Yang et al. Solar Physics, 283, 
369, 2013). This method requires that the magnetic flux to be balanced on all the side 
boundaries of the considered volume. In this paper, we propose a scheme to obtain 
the vector potentials at the boundaries to remove the above restriction. We also used a 
theoretical model (Low and Lou, Astrophys. J. 352, 343, 1990) to test our scheme. 

Subject headings: Magnetic helicity 
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1. Introduction 



Magnetic helicity is a key geom etrical parame ter to describe the structure and evolution 



of solar coronal magnetic fields ( e.g. iBerger 
determined as 

where A is the vector potential for the 



19991). Magnetic helicity in a volume V can be 



Jv 



A ■ BdV, 



(1) 



magnetic field B in this volume. Magnetic helicity is 



conserved in an ideal magneto-plasma (|Woltjei 



1958h . As long as the overall magnetic Reynolds 



number is large, however, it is still approx imately conserved, even in the course of relatively 



slow magnetic reconnection (|Berger 



19841) • The concept of magnetic helicity has successfully 



been applied to characterize solar coronal processes, fo r a recent review about mod eling and 
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tOO% . Despite of 



observations of photo spheric magnetic helicity see, e.g., 

its important role in the dynamical evolution of solar plasmas, so far only a few attempts have 

been made to estimate the helicity of coronal magnetic fields based on observations and numerical 

simulations (see, e.g., Thalmann, Inhester, and Wiegelmann, 2011; Rudenko and Myshyakov, 

2011). 



Yang et al. 



(|2013h developed a method for an efficient calculation of the rela tive magnetic 
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helicit y in finite 3D volume which already was applied to a simulated flaring AR 
(|201l|) . This method requires the magnetic flux to be balanced on all the side boundaries of the 
considered volume. In this paper, a scheme to remove the restriction has been proposed. In Sec. l2l 
we describe the restriction of vector potential in the former paper. In Sec. |3l we present the details 
of the new scheme to calculate the vector potentials on the six boundaries. In sec. |4l we use the 
theoretical model to check our scheme. The summary and some discussions are given in Sec. |5l 



2. The former definition of Ap and A at the boundaries 

Let us define a finite tliree-dimensional (3-D) volume ("box") in Cartesian coordinates with a 
magnetic field B(jic, J, z) given in this volume. Let the volume be bounded by x = [0,lx],y = [0, /j], 
andz = [0,Z,]. 

First one has to provide the values of Ap and A on all six boundaries {x = 0,lx',y = 0, ly; z = 
0, 1 J. To take the bottom boundary (z = 0) for example, we define a new scalar function (p{x,y) 
that determines the vector potential Ap of the potential magnetic field P on this boundary as 

follows: 

dip dip I 

According to the definition of the vector potential, the scalar function ip{x, y) should satisfy the 
Poisson equation: 

^ip{x,y) = B,{x,y,z = Q). (3) 

The value of d(f>/dn on the four sides of the plane z = is set to zero in Equation (|3). According 
to Eq.([2l), Apjc and Ap,, will vanish aty = 0, ly and at x = 0, l^, respectively, on the z = plane. 
Thus, the corresponding magnetic flux at the boundary should also vanish because of Ampere's 
law. The values of Ap on the other five boundaries could be obtained in a similar way. For the 
vector potential A at all boundaries the same values are taken as for Ap. When the magnetic fluxes 
at the six boundaries are not zero, we should calculate the value of vector potentials at the twelve 
edges of the three-dimensional (3-D) volume to provide the Neumann boundary for the Poisson 
Equation at each side boundary. In next section, we will introduce a scheme to calculate the vector 
potentials at the twelve edges. 




Fig. 1. — Magnetic flux O,- (/ = 1, ..., 6) at the six boundaries and the integrals a, (/ = 1, ..., 12) of 
J Ap • dl at the twelve edges. 

3. new scheme to obtain A^ and A at the boundaries 

For the Ap, we define the magnetic flux O, (/ = 1, ..., 6) respectively at each side boundary 
(z = 0; z = Iz', X = 0; x = I/, y = 0; y = ly). The integrals of jAp ■ dl at the twelve edges are 
defined as a, (?' = 1, ..., 12). The twelve integrals and the corresponding directions are represented 
in Fig. [B 

According to the Ampere's law, the integral value a, satisfy the following linear equations 



TX = B, 



(4) 



where B = (01,02,03,04,05,05)^, X = {ai,a2,a-i, ...,^12)^ and T is a matrix of 6x12 , which is 
equal 



111100000000 

000011110000 

000-1000-11 00 1 

0-1000-1000 110 

-10 0-10 0-1-10 

0-10 0-10 0-1-1 
One can check that the six rows-vector in this matrix is not linear independent because the 

magnetic field is divergence free and the sum of O, at the six boundaries is zero. Moreover, the 

unknown twelve a, are not unique just under the restriction of above six conditions. Hence, we 

need to construct twelve independent conditions to obtain the unique solution for a,. We define 

the new matrix T as follows 

111100000000 

000011110000 

0-10 0-11 001 

0-10 0-10 110 



(5) 



■10 0-1000 



-1-10 



1 







1 



-1000000000 







-100000000 

00 10-10000000 

000 10-1000000 

10-100010 

0-10 

1 0-1000 

One can check that the determinant of T is not zero. According to Cramer rule, the unique solution 
is existent for the new linear equation 

tX = B, (7) 





1 





(6) 



where X = (ai,a2, a^, ..., auY and B = (<I>i, <I>2, <D3, <D4, O5, 0, 0, 0, 0, 0, 0, 0)^. Then we can obtain 
the integrals of jAp • dl at the twelve edges. The corresponding vector potential at the twelve 
edges could be obtained by using the following equation: 



Apx ia,) = — sinCTTx/L.v), / = 1, 3, 5, 7 

Apy {ad = — - sininylLy), i = 2, 4, 6, 8 (8) 

Apz {ad = ^ sin(;rz/L,), / = 9, 10, 1 1, 12 

Note that Ap at the ends of every edge both are zero according to the above equation. That 
is the requirement of Eq. (|2]). Then we resolve the Poisson equations to obtain Ap at the six 
boundaries. For the vector potential A at all boundaries the same values are taken as for Ap. Then 
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(120 13h to calculate 



we can follow the method of Sec. 2.2 and 2.3 of the former paper 
the relative magnetic helicity in this volume. 

4. Testing the scheme 

For testing the new scheme to obtain the vector potentials at the boundaries, we use the 
axisymmetric nonlinear force-free fields of Low and Lou (1990). We used the model labeled Pi 1 
with / = 0.3 and O = n/l in the notation of their paper. We calculated the magnetic field on a 
uniform grid of 64 x 64 x 64. The pixel size in the calculation is assumed to be 1. 

We calculate the magnetic fluxes <I>o at the six boundaries and substitute it to the Eq. (|7]) to 
obtain the integral a, at the twelve edges of the 3D volume. Then we substitute a, into Eq. ^ 
respectively to get the boundary value for resolving the Poisson equation in Eq. © at the six 
boundaries. After we attain Ap at the six boundaries, we could also calculate the magnetic flux O 
according to the relation between the vector potential and the magnetic field: B ■ n = V x Ap • n. 
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Table. [1] represents the final result after we apply the the above scheme. It can be found that 
the calculated magnetic fluxes at the six boundaries by using our scheme respectively coincide 
well with the original value from the theoretical model. Note that the total magnetic flux of the 
theoretical model is not exact zero. However, it is required that the total magnetic flux is exact 
zero when resolving the linear equation Eq. (|7}, which cause the total magnetic fluxes of ^ Ap ■ dl 
and O are diff"erent with that of Oq. On the other hand, the numerical errors when resolving 
the Poisson equation are also unavoidable, which will also introduce the difl"erence for the total 
magnetic flux as well. 

5. Summary 

In this paper, we propose a new scheme t o calculate the vec tor potential at the boundaries 



to remove the restrictions in the former paper lYang et al. 



(|2013h . In principle, now we can 



calculate the relative magnetic helicity of any magnetic field structure in Cartesian coordinates. In 
the observations, we could use force-free extrapolation method to obtain the three-dimensional 
magnetic structure to analyze the evolution of relative magnetic helicity. On the other hand, we 
can also use a seq uence of magnetograms to estimate the accumulated magnetic helicity in the 



solar corona (Ref. 



Demoulin and Pariat 



I2009h . It will be very interesting to compare the two types 
of accumulat ed magnetic helic ity and analyze the correlation between magnetic helicity and solar 
eruption (e.g. Jing et al.ii2012,) . In the simulations, we could also calculate the relative magnetic 
helicity directly based the known magnetic field structure to understand how the magnetic helicity 
plays an important role in solar reconnection and dynamos. 

This study is supported by grants 10733020, 10921303,41174153,11173033 11178016 and 
11103038 of National Natural Science Foundation of China, 201 1CB8 11400 of National Basic 
Research Program of China, a sandwich-PhD grant of the Max-Planck Society and the Max-Planck 



Society Interinstitutional Research Initiative Turbulent transport and ion heating, reconnection 
and electron acceleration in solar and fusion plasmas of Project No. MIF-IF-A-AERO8047.The 
authors also like to thank the Supercomputing Center of Chinese Academy of Sciences (SCCAS) 
for the allocation of computing time. 
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Table 1 : Testing using the new scheme to a theoretical model. 



side boundary 


a 
^0 


^Ap . dl' 


c 


z = 


-3615.81 


-3615.81 


-3487.1068 


z-k 


1461.13 


1461.13 


1490.2739 


x = 


-1471.95 


-1471.95 


-1563.7657 


x = h 


-1471.95 


-1471.96 


-1564.9280 


y = 


4006.42 


4006.42 


4087.8426 


y-iy 


1068.27 


1092.17 


1066.6782 


Total flux 


-23.89 


0.00012 


28.99 



"Magnetic flux at side boundary of theoretical model. 

''The integral of Ap ■ dl for each side boundary. 

^Magnetic flux from the solution in Sec. |2]at each side boundary. 
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